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Abstract 



We construct the higher twist structure functions that describe inclusive b 
hadron decays in regions of phase space where the hadronic decay products 
carry high energy but have low invariant mass. We show that, for B meson 
decays, there are four new non-vanishing matrix elements of non-local oper- 
ators. We show that to subleading twist these decays are parametrized in 
terms of four functions. We compute the tree- level matching for a general 
heavy-to-light current and apply it to B — ► X s j. Using a simple model for 
these functions we estimate the subleading twist contributions to this decay. 
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I. INTRODUCTION 



Inclusive decays of B mesons can be analyzed via an operator product expansion (OPE) 
HI as a power series in AQ CD /m b . At leading order the OPE reproduces the parton model 
results, while non-perturbative physics is parametrized by matrix elements of higher dimen- 
sional operators and suppressed by powers of the heavy quark mass. 

In practice, higher dimensional operators in the OPE are only suppressed for sufficiently 
inclusive observables. In the resonance regime, where the invariant mass of the final hadronic 
state is restricted to be <Aq CD , the decay is no longer inclusive and it is not surprising 
that the OPE completely breaks down. However, the OPE includes terms suppressed by 
-^xAqcd/^x (where Ex and mx are the energy and invariant mass of the final hadronic 
state) which are suppressed over most of phase space, but are 0(1) in the region of high 
energy, low invariant mass hadronic states, 

E x ~ m b , m\ ~ AQ CD m 6 > A^ CD . (1.1) 

In this "shape function region" the OPE breaks down, even though it is far from the reso- 
nance regime. Thus an inclusive description is still appropriate, but the expansion parameter 
in the OPE has to be modified. It has been shown that the most singular terms in the OPE 
may be resummed into a non-local operator [[| 

O (u) = 7iJ(uj + in ■ D)h v , (1.2) 

where is a light-like vector in the direction of the final hadrons, h v is a HQET heavy 
quark field and we denote variables normalized to m& by a hat: = D^/rrib- The matrix 
element of this operator in a B meson is the light-cone structure function of the meson, 

f(tu) = ^(B\O (u)\B). (1.3) 

The rate in the shape function region is determined by f(oo). However, since it is a non- 
perturbative function, f(u) cannot be calculated analytically, and the rate in the shape 
function region is model-dependent even at leading order in AQ CD /m b . 

Unfortunately, for B — > X u iv decay (and to a lesser extent B — > X^), the experimental 
cuts which must be imposed on the phase space to eliminate large backgrounds from b — > 
c decay typically put the decay into the shape function region, introducing large model- 
dependence in the predicted rate. For semileptonic b — > u decay, this is the case for cuts 
on either the charged lepton energy or the hadronic invariant mass ||, and this model 
dependence is the major theoretical stumbling block to a precise determination of the CKM 
matrix element \V u b\ from inclusive decays. 

However, since the distribution function (|1.3|) determines the shape of the photon spec- 
trum in B — ► X s 7 as well as the charged lepton or hadronic invariant mass spectrum in 
B — > X u tv, it was suggested a number of years ago J|] that f(y) could be measured in 
B — ► X s 7, and then used to extract \V u b\ from semileptonic decayQ The perturbative cor- 



* Another solution is to consider a kinematic cut which does not put the final state into the shape 
function region, such as a cut on the lepton invariant mass ||. 
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rections to the relation between these processes have been intensely studied in recent years 

In addition to the radiative corrections, however, there are non-perturbative corrections 
to the relation between the photon spectrum in B — > X s 7 and the lepton or hadronic 
invariant mass spectra in B — > X u iv. The light-cone distribution function ( |1.3| ) only resums 
the most singular terms of the OPE in the shape function region. There are corrections to 
this from less singular terms suppressed by Aqcd/^t-6, analogous to higher-twist corrections 
to deep inelastic scattering J7j . These corrections are important for a precision measurement 
°f \Vub\ y i a this method, but have not yet been studied. 

In this paper we discuss the subleading corrections to heavy-light B decay in the shape 
function region by performing a twist expansion rather than the usual OPE in terms of local 
operators^. At leading order in the twist expansion we reproduce the known results, while 
at subleading order we find four new non-local operators relevant for B decays. We compute 
the tree-level matching for a general heavy-light current and apply it to B — > X s 7, reserving 
a discussion of B — > X u tv decay for a future work 0. We use a simple model for these 
functions to estimate the subleading twist contribution to B — ► X^. 



II. KINEMATICS 

We consider a general heavy to light transition (radiative or semileptonic) , proceeding 
via the current 

j(x)=q(x)Tb(x), (2.1) 

where q(x) is a massless quark field and T is an arbitrary Dirac matrix. The decay rate is 
related to the imaginary part of the T-product of two heavy-light currents: 

dT 21m(B\T\B), (2.2) 

where 

T(q) = f*xT[j(0)ji(x)]e i «', (2.3) 

and q is the momentum transfer. 

The kinematics for this process are shown in Fig. |l[ As usual, the heavy quark momentum 
is split into a large and a residual piece, = + k^. In the shape function region the 
final hadronic state has large energy but small invariant mass, and so its momentum lies 
close to the light-cone. We introduce a light-like vector n M to define the expansion of the 



t As discussed in [||, there are actually two stages of matching: at \i = QCD is matched 
onto an intermediate theory with collinear and soft degrees of freedom, while at a lower scale the 
non-local OPE is performed. Since we are not concerned with summing Sudakov logarithms in this 
paper, we may neglect the intermediate theory. 
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p b = m b v + k 



FIG. 1. Kinematics for a general heavy to light transition. 

jet about the light-cone, and the velocity of the heavy quark then defines a second light-like 
vector n = 2v — n. These vectors satisfy 

n 2 = = n 2 , v ■ n = v ■ n = 1, n ■ n = 2. (2-4) 

In the frame in which the b quark is at rest and the emitted hadrons move in the +z direction, 
these vectors are given by n M = (1,0,0,1), = (1,0,0,-1) and t> M = (1,0,0,0). We can 
decompose the metric according to 

= - (n»n u + n u n») + , (2.5) 

which defines g^ . 

In the shape function region, the momenta scale as 

[m b v — q) ■ n = m b — q ■ n ~ 0{m b ), 

(m b v - q) ■ n = m b - q ■ n ~ 0(A QCD ), (2.6) 
fc" ~ 0(A Q cd). 

It is therefore convenient to split the momentum q 11 into large and small components, 

(f = qm + r ) (2.7) 

where 

Q» = \(m b n^ + q-nn fl ), F = -\{m h - q ■ n)n^ (2.8) 

are 0(m b ) and O(Aqcd), respectively. The momentum of the light hadronic decay products 
is 

p v = m hV » - Q» - p + W = \{m b - q ■ n)n^ + \(m b - q ■ n)n^ + W , (2.9) 



and so in the shape function region ( |2.6| ) we have p = C^Aqcd^) *C rn 2 . 



III. MATCHING 
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A. Leading Order 

The expansion in the shape function region differs from the usual l/m^ expansion because 
of the additional small parameter m& — q ■ n; in the usual OPE terms of order 



m& — q ■ n 



are treated as subleading, whereas they are 0(1) in the shape function region. Thus, instead 
of resumming terms of this form to all orders, we perform an OPE in powers of Aqcd/^;), 
but using the scaling (|2.6|) . This is analogous to the twist expansion in DIS. 

Expanding the light-quark propagator shown in Fig. [I] in powers of AQ CD /mb gives 

m h j + 1 1/1 ^ /^Aqcd 



O . (3.1) 



{m h v + k - Q - l + ie) 2m b -q-n + k- n + ie V m b 

Since k ■ n and m b — q ■ n are of the same order, this term cannot be expanded in powers of 
k ■ n, and so cannot be matched onto a finite set of local operators. 
Instead, consider the set of operators 

6 {t) = h v {0)E{0,t)h v {t), (3.2) 

where we use the shorthand notation 

V>(t)=V>(— ), (3.3) 

where t is dimensionless, to denote fields on the light-cone defined by n M . The path-ordered 
exponential 

E(tx, t 2 ) = P exp (-i ^ n ■ A(t') dt'^j (3.4) 

is required to make the operator Oo(t) gauge invariant. The operators Oq(uj) are defined in 
terms of Oo(t) by the linear combination 

1 r°° 

O (uj) = — dt e~ luJt O {t) = h v 5(u + in ■ D)h v , (3.5) 

and have the required form for the imaginary part of the leading term in the heavy quark 
expansion Q3.1| ). The matrix elements of Oq(u) define the light-cone distribution function 
of the b quark in a B meson 

/(o;) = — *— (B\h v 8(cu + in ■ D)h v \B) . (3.6) 
Expanding f(uj) in powers of in ■ D gives the series of increasingly singular terms 



Al r,/, N Pi - w , 
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/H = «M- 5 V(«)-£H-H + - (3.7) 



where 



■(B\h v (iD a )(iDp)h v \B) = -(g a/3 - v a vp)\ 1 



2m b 

-^—{B\{iD a ){iD^)(iDp)\B) = ^{g a p - v a V/j)v M pi. 



(3.8) 
(3.9) 



For a general h hadron decay there is an additional parity-odd operator at leading twist 

P«(t) = h v (0)j a ^E(0,t)h v (t) , (3.10) 

whith Pq(u) defined analogously to ( |3.5|) . This operator is not relevant for B meson decays 
since its matrix element vanishes, but it gives a spin dependent contribution to A& decay. A 
general Dirac structure between heavy quark fields may be expressed in terms of these four 
independent matrices via the projection formula 



ip + Tr(P + r) - ^Tr(sT) , 



where 



+ • 



(3.11) 



(3.12) 



andP+ = ±(l + ?4). 

Since the OPE is performed over a a continuously infinite set of operators labeled by u, 
the heavy quark expansion in the shape function region is, like in DIS, a convolution over 
a single parameter which may be interpreted as the light-cone momentum fraction of the 
heavy quark: 



ImT(g) 



2rrih J- 



duo 



C (lu)O (u) + C? (u)P 0>a (u) + O 



V)(T>\ 



(3.13) 



where the Cj's are perturbatively calculable short distance coefficients. The tree level match- 



ing conditions are easily obtained from (3T) and ( 3.11 



C (v,q,uj) 



7T„ 



THP+r^rWl-n-g-w) 



7T, 



Tr(s Q r^r) 5(l-n-q-u) 



(3.14) 
(3.15) 



B. Subleading Order 



Expanding ( |3.1|) to subleading order in 1/rrib, we will in general match onto non-local 
objects of the form |7| 

O^-^ih, ...,tn) = K^TiiD^ih) . ..iD^(t n ^)]h v (t n ) , (3.16) 

where T = 1 or 7^75 and 
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D„{t) = d„ - igA a Jt)T a 



(3.17) 



is the usual covariant derivative, acting at a light-cone coordinate. 

The most general case involves operators in which every field and derivative is evaluated 
at a different light-cone coordinate. However for heavy quark decays to subleading order, we 
find that at maximum only two light-cone coordinates enter. The complete set of operators 
required is: 

or(*i,t 2 ) 

P^itiM) 
6f(ti,* 2 ) 

Pla(tl,h) 



h v (Q) (-iD' 1 (Q)E(0,t) + £7(0, t)iD"(t)J h v (t) 
h v (0) (-W^(0)E(0, t) + E(0, t)iD»{t)) lal5 h v (t) 
ih v (0) (-iD»(0)E(0,t) - E(0,t)iDi*(t)) h v (t) 
ih v (0) (-W>*(0)E(0,t) - E{0,t)iD"(t)) WsKW 
h v (0)E(0, *0 {iDtih), iDHh)} E(t u t 2 )h v (t 2 ) 

^(0)S(0,t0{^(t0,^I(^)}^(^l ) *2)7a75^fe) 

gh v (0)E(0, tOGf (t0^(ii, t 2 )h v (t 2 ) 



gh v (0)E(0, t x )Gf (h)E(ti, t 2 ) lalh h v (t 2 ) , 

where Dj_ = g^D v and gG^ = i[(iD±), (iDjJ\ is the gluon field strength. 
The Fourier transformed operators Oi are defined as 



<9 2 » 



Ot v {u x ,u 2 ) 



— / dt e- mt O? (t) = h v {«£>", 5(in ■ D + u) } h v 

iD<*,6(in- D +u. 



oo 
oo 



2 

dh dt 2 e~ iuJltl 



ih, 



1 

2^ 



h v S(in ■ D + io 2 ) {iD^_, iD v ± } 5{in ■ D + u)i)h v 

1 \ poo _ 

— / dt.dhe-^e-^-^O^it^h) 

III J J-oo 



= gh v 5(in ■ D + uj 2 )G^ 5(in ■ D + u)\)h v 
Similarly the Fourier transforms of the Pj's are 

P£Jw) = h v {iD^ 5(in ■ D + oo)) lalb h v 



P 2 , a (uj) = ih v iD**, 6 (in ■ D + 10) 7^75/^ 

h v 5(in ■ D + uo 2 ) {iD^_, iD 1 ^} 5(in ■ D + uoi)^j a ^h v 
P£a(ui, uj 2 ) = gh v 8(in ■ D + u 2 )G^8(in ■ D + a;07a75^- 



(3.18) 



(3.19) 



(3.20) 



The Feynman rules for the operators Oq — O4 in n ■ A = gauge are shown in Fig. ^. 

Finally, at subleading order there are also contributions from the time-ordered products 
of Oq(lv) with the subleading terms in the HQET Lagrangian, 
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S(u> + n • k) 




2^5^ + n-k) 






a, a 

-gT a g' ia [5(u + k ■ n) + 

8{u + n ■ (fc+p))] 

Of (wi,w 2 ) 
<5<j 



a, a 

-iffTV"* [5(lu + k ■ n) - 

<5(w + n ■ (k + p))] 

Of (wi.wg) 




-, 9 T« [(p + 2k)lrf" + (p + 2*)isS°] x 
<5(a,'i + n • fc)d(aj 2 + n • (p + fc)) 






Pi a " a ft 6 P2 



Pi a ,3, 6 P2 

*g 2 {gTgf-gfsT)\T a ,T b ) 



x 5(uj 2 + n ■ (k + pi +p 2 )) S{ix>i + n ■ k)S (ui 2 + n- (fc+pi +P2)) <5(^i + n • fc)i5 (w 3 + n • (fc + pi +P2)) 

FIG. 2. Feynman rules for non-local operators Oo — O4 in n • A = gauge. The Feynman rules 
for Pq — P4 are identical except for the Dirac structure. 
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C 1/m (y) = h v (y)(iD) 2 K(y) + ^h v (y)a^h v (y) . 



This yields another two operators 



OrW) > / d 4 y i- Jdte-^T (h v (0)h v (t)O 1/m (y) 



(3.21) 



(3.22) 



At subleading order the nonlocal OPE in (|3.13|) is 

-2m 6 ImT(g) = fdoo (C (v, q, u)O q (oo) + C^ (v, q, u)P 0>a (. 

Am b i=\o J 



(00 



■ — V [duJxd^ (C^(v,q,u}x, 002)0^(^,002) (3.23) 



2m b i=3>4 

+Cgf(v, q, L0 X , UJ 2 )Pi,a,^(Ul, CO2) 

1 r 

— doo [C T (v, q, uj)O t (uj) + C* >T (v, q, oj)P T>a (u) 



2mb 
+0 
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The matching at subleading order onto the operators ( |3.18|) is performed by computing 
the zero, one and two gluon matrix elements of (|2.3| ) in full QCD and comparing this to the 
operators in ( |3.19| ) and ( |3.20| ). Note that this includes terms from the expansion of the b 
quark field, 

6 =( i+ S + -K p- 24 » 

At tree level, we find 

C?(v, q, u) = ^Tr ({P + , -f}^) 5(1 - n ■ q - u) 
C 2 >, q, u) = (-z)yTr ([P+, -f]Ti0) 5(1 - n ■ q - u) 



4 

Cn«,g,w 1)Wa ) = ~gT — Tr (p + r^r) (3.25) 

2 m b — n ■ q v ' 

'5(1 — n ■ q — oji) — 5(1 — n ■ q — uj 2 )\ 



x 



UJ\ — u 2 



Cr(., ?)Wl , W2 ) = 4 -TrfPT^-Or) 

2 m b — n ■ q v 1 

( 5il — n ■ q — uji) — 5(1 — n ■ q — uj 2 ) 



V LU 1 —U 2 J 

C T (n -q,u) = C (n ■ q,u) , 
and for the corresponding spin dependent operators 

Cgf («, q, u>) = -^Tr ({ a «, YWT) 5(l-n-q-u) 
C${v, q, w) = i ^Tr 5(1 - n • £ - w) 

C^f>,<Z,^ 2 ) = — Tr (s«TiT) (3.26) 

2 777ft — n ■ q K ' 

( 5(1 — n ■ q — uji) — 5(1 — n • q — 0*2) \ 

\ U1-U2 ) 

CgT(n ■ g,^,cu 2 ) = (-<)£ — Tr ( S ^(-<)r) 

2 mft — n ■ q K J 

5(1 — n ■ q — o^) — 5(1 — 77 • <J — 0^2) \ 



x 



\ CJi — u 2 

C? T (n -q,u) = C? (n ■ q,u) . 



IV. MATRIX ELEMENTS 



Matrix elements of the subleading operators ( |3.19|) , ( |3.20| ) and (|3.22| ) give rise to new, 



subleading structure functions. Writing the most general Ansatz consistent with the sym- 
metries and the equation of motion (iv ■ D) h = 0, we find that only the following matrix 
elements are non-vanishing 
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(B(v)\0^)\B(v] 
(B(v)\0^(u l7 u 2 )\B(v) 
(B(v)\P 2 %(u)\B(v] 

(B{v)\Pj^{oh^)\B(v) 
(B(v)\0 T (u)\B(v] 



2m B - n") 

2m B g 2 {uj 1 ,u 2 )g± / 
2m B /ii(w)e^ a 

2m B h 2 {ux, uj 2 )e paa p g^g v °v p 
2m B t(u) , 



(4.1) 



where we define 



(4.2) 



and e 0123 = 1. 

The matrix element of 0^ u (ui, u 2 ) between B mesons vanishes since no antisymmetric, 
parity even object can be constructed which is perpendicular to both v and n. Similarly, 
the matrix element of P^ v a {ui,u 2 ) vanishes, since there is no parity odd, symmetric object 
perpendicular to v and n. Due to the equations of motion, the matrix element of 2 {u) 



must be proportional to {v^ 

(B{v)\0^{uj)B{v)) = - n") . 
Contracting with n M we find 

a= (B(v)\n ■ 2 (u)\B(v)) = (B(v)\ih v [in ■ D, 8(in ■ D + uj)]h v \B(v)) = 



(4.3) 



(4.4) 



(where the last equality is due to the delta function) and so the matrix element vanishes. 
The matrix element of P\ a (u}) vanishes since all of its moments with respect to uj vanish. 
Finally, the matrix element of Pr >Q (o;) vanishes due to parity. 

There is additional information on the remaining new functions. Starting with g%(uj), we 
find 

2m Bgi (uj) = n^{B(v)\0^(uj) \B{v)) = (B(v)\h v {in ■ D, 5 (in ■ D + cu)} h v \B(v)) 

= -2(m b uj)(B(v)\h v 5(iD ■ n + u)h v \B(v)) = -Am B (m h uj) f(u) . (4.5) 

Thus, gi(u) is determined by the leading order structure function, gi(u) = —2{m b uj) f(u). 
Some information can also be obtained on moments of the function g 2 (uJi,uj 2 ): 



(m,n) 
92 



(4.6) 



2m 



(B{v)\h v {in ■ D) m (iD ± ) 2 {in ■ D) n h v \B(v)) 



B 



leading to 



g .f» 

92 



2\l 
~3~' 

„(0,n) 
92 



(4.7) 



where the last equality arises because of the constraints from the equations of motion 
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(B\h v {iD a ){iD ul ) . . . {iD Un ){iDp)h v \B) = {g af} - v a v p )A Vl .. Mn . (4.8) 

We can also obtain information on the parity odd operators. The function hi(u) is a genuine 
new non-perturbative function which introduces spin dependent effects. The first three 
moments of this function are given by 

[duj hi (oj) = ; I duo u hUu) = -; fdu^hiiuj) = (4.9) 

J J nib J m t 

where 

-^—{B\h v {iD a ){iD p )s x h v \B) = \ie uapx v v \ 2 (4.10) 

Zui b Z 

1 1 

{B\h v {iD a )(iD ll )(iDp)sxh v \B) = -ie ua pxv u v^p 2 . ( 4 -H) 

The function h 2 (wi,uj 2 ) also introduces spin dependent effects with the first few moments 
of the function given by 

/4°' 0) = A 2 ; ^ 0) = = ^°' 1) . (4.12) 

Finally we have to consider the function t(uo). Since l/m^ terms are absent in the total rate, 
we have 



diut(Lu) = 0. (4.13) 

Furthermore, the first moment of t(u>) is related to Ai and A2 

2m B J dcu to t(cu) = -i j d A x (B\T {[h v (in ■ D)h v ](0) 1/m (x)} \B) (4.14) 

Ai + 3A 2 



2m 



b- 



while the second moment introduces one new parameter 



duuj 2 t(to) = , (4.15) 



where 3r = —2(71 + 2>T 2 ) and % are the parameters used in ]TT],|T2 



Combining these results with the known leading twist contribution ( |3.7|) leads to 



/nM = — <S'M + + (4.16) 

rrib zmi 



h 2 {uj l ,uj 2 ) = X 2 S(uj 1 )5(uj2) H 

m b zmf 
11 



where we have used the relations 

re = -2*'({); = -35"(0 , 

which are true when integrated against a function which is non-singular as £ — ► 0. 



V. APPLICATION TO 5 —> X s j 

The decay i? — > X s 7 is described by the effective Hamiltonian 

H c f[ ■■ 



.^ Vtb v; s cMo 7 + . 



where 



On 



32tt ; 



and the dots denote additional operators which we shall neglect for the purposes 
discussion. This effective Hamiltonian leads to the Dirac structure 

r = -iAa^e,q v (l + 75) with A = ^V tb V*C 7 (^)^m b . 

The kinematics of this decay are particularly simple since q 2 = 0, 

q^ = E^-nT = x—rf 

and so the large and small kinematical factors defined in section || are 
m b — n ■ q = m b , m b — n ■ q = m b (l — x) . 
Computing the matching coefficients, we find for the rate 

To {/(I -x)- 2(1 - x)f(l -x) + -U( 



dT 

dx 



[l-x) + — hAl-x) 



m b 



mi 



G 2 (l - x) - H 2 (l - x) 



where 



and 



G 2 F a\V ts V t t\ 2 \C 7 (m b )\' 
32vr 4 



-m 



b i 



G 2 {l-x) 
Ho(l-x) 



dui duj 2 g 2 (uj 1 ,uj 2 ) 
duji duj 2 h 2 {ui,u 2 ) 



5(1 — x — uji) — 5(1 — x — uj 2 ) 
5(1 — x — uji) — 5(1 — x — uj 2 ) 

Ll>i - UJ 2 



(5.1) 

(5.2) 
of this 

(5.3) 

(5.4) 
(5.5) 



(5.6) 

(5.7) 

(5.8) 
(5.9) 
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This result can be compared with previous work [ 12 1 , where the 1 / m& expansion has been 
discussed. This expansion is recovered by performing a moment expansion of (|5.6|) . Using 
( CT ) we find 



2Ai 



6'(1 



3 

-X 2 6'(l-x) 



x) 



leading to 

dT 

dx 



Tr 



6(1- 



x) 



G 2 (l-x) 
H 2 (l - x) 

_ Ai + 3A 2 
Ai 2pi - 3p 2 + Tx + 3T 2 



S'(l-x) 



(5.10) 



(5.11) 



6m? 



+ 



6m? 



5"{l-x) 



Pi 



18m? 



6"'(1 



x 



in agreement with [F2|. From the matching coefficients ( |3.25|) we see that the Wilson co- 
efficients Ct is identical to the leading order coefficient. This implies that for any current 
mediating a heavy-to- light decay the subleading structure functions t(u) always arises in 
the same linear combination with the leading order function f(uj). Thus, we can always 
combine the functions fiui) and t(u) into a new universal function which is defined by 



F(u) = f(u) + - — t(c^) 



2m h 



This new structure function has the moments 

Xi + 3A 2 



duF(u) 



1; 



duo ijj F{u) 



2m 2 b 



dujuj F(uj) 



Ai 



3m^ 



+ 



2mf 



(5.12) 



(5.13) 



There are three new subleading twist structure functions; the spin independent function 
G 2 (u), as well as hi(oj) and H 2 (uj), which are sensitive to the heavy quark spin. Thus there 
are in total four functions parametrizing the heavy-to-light decays to subleading twist, since 
we may replace 



uifiuj) = oj F{bS) + . . . 
where the dots denote higher twist terms. 



(5.14) 



A. A Simple Model 



To get some insight into the size of the effect of these new functions we will use a simple 
model for these functions which incorporates the information we have on their moments. 
The following model has been proposed H for the leading twist distribution function 



, / x 32m fe 

/modM = ~v [1 



7T 



yfexp 



7T 



e(i-v) 



y 



(5.15) 



where A is the only free parameter. This model therefore assumes a simple correlation 
between all higher moments of f(oo). 
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All of the new subleading twist functions have vanishing zeroth moments. In order to 
construct a model for these subleading functions we use the derivative of the leading twist 
function. Normalizing this derivative to match the known first moments of the subleading 
functions, we obtain 



•fmod(k') — /mod^) 
^2 mod (^0 



A + 3A 2 



2m\ 



/modM 



2Ai , , v 



(5.16) 



For the decay B — > X s ^y the two functions hi(u) and H2{uj) enter in the combination 

A h {uo) = h l {uj) + H 2 {uo)/m b . (5.17) 
The zeroth and first moments vanish, while it has a non-vanishing second moment 



du uj 2 A h (uj) = . 

mi 



(5.18) 



In our simple approach to modeling the subleading functions, we would obtain A fe mod (u;) = 
0, since the first moments of m b hi{uo) and H 2 (u!) coincide. To use the information on the 
second moment ( |5.18| ), we instead model A h (u) by the second derivative of / mo d : 



A^modO^O 



P2 

2m\ 



•/'in 



mod 



to 



This leads finally to our model for the differential decay spectrum of the decay B 

cT'mod 



dx 



To 



(2X - l)F mo d(l ~X) + |^/mod(l ~ x) + £jpfLi 




mod I 



X) 



3m b 
_ Ax + 3A 2 
2m\ 



/mod(l X ) 



+ 



1 — X 

2Ai 
3mt 



(5.19) 

x sT . 

(5.20) 



/mod(l X ) 



2m^ mod ' 



l-x) 



We can now use this spectrum to analyze the effect of the subleading twist contributions 
to the partially integrated decay rate 



He") 



i 



— dx 

dx ' 



m B 



and xq 



2E^ 
m b 



(5.21) 



with E® being a lower cut on the photon energy. The effects of the subleading shape 
functions are shown in Fig. |3|, in which we plot the ratio of the partially integrated rate 
with and without the subleading twist contributions as a function of the photon energy cut 
for various values of the parameters A and p 2 . With this simple model, the curves on this 
plot should only be taken as an estimate of the size of the corrections in different kinematic 
regions. Recall that we expect that for a photon cut E® <ms/2- Aqcd ~ 2.1 GeV the usual 
OPE should hold, while for m B /2 - A QCD ~ 2.1 GeV<E°<m B /2 - A 2 QCD /m B ~ 2.5 GeV 
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T ( E y)leading 0-9 
twist 

0.8 
0.7 
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Ey(GeV) 

FIG. 3. Partially integrated Rate ( |5.21| ), normalized to the leading twist result, using the 
simple model given in ( |5.20|) . The solid, short-dashed and long-dashed lines correspond to A = 570 
MeV, 470 MeV and 370 MeV, respectively. The lines which rise at the endpoint correspond to 
P2 = (500 MeV) 3 , while those that go down correspond to p2 = —(500 MeV) 3 . The values 
of Ai has been chosen to reproduce the second moment of the leading order structure function, 
Ai = -0.53 A 2 . 

the twist expansion presented in this paper is appropriate, with subleading twist corrections 
naively of order AQ CD /m b ~ 10%. For E®>rriB/2 — AQ CD /mB ~ 2.5 GeV we are in the 
resonance region, and the twist expansion is expected to break down. From the figure we see 
that these expectations are borne out in our simple model: below ~ 2.1 GeV the corrections 
to the leading twist are small, since they are only corrections of order AQ CD /m 2 to the 
leading rate. Below E® ~ 2.4 GeV, the corrections are less than or of order ±20%, while 
the twist expansion starts to break down above this cutoff. 

VI. CONCLUSIONS 

It has been known for some time that in inclusive heavy hadron decays the naive short 
distance expansion must be replaced by a twist expansion if the phase space is restricted 
to a region of large energy, low invariant mass final hadronic states. The leading term, 
parametrized by the the light-cone distribution function of the heavy quark in the hadron, 
is well investigated by now, but subleading terms of this expansion have not been previous 
studied. In the present paper we have identified the non-local operators appearing at sub- 
leading order in the twist expansion. The tree level matching to these operators has been 
computed for a general bottom hadron decay and the matrix elements of the subleading 
operators have been parametrized for B meson decay. 

We found that for any inclusive B meson decay four independent subleading distribution 
functions are needed. We worked out the case for B — > Xg'j in detail. Using a simple 
model for the leading and the subleading distribution functions we studied the effects of the 
subleading terms on the photon energy spectrum. We found that they had the expected 
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behavior: in the region where the local OPE is appropriate, these corrections were negligible, 
whereas in the region where the twist expansion was appropriate, they were of order 10-20%, 
depending on the parameters of the model. 

Since the leading distribution function is not known, much less to 10-20% accuracy, these 
results are of limited utility for B — ► Xgj decays (although they do indicate the region where 
the twist expansion breaks down). However, there are certain relations between the charged 
lepton energy spectrum in B — > X u lv and the photon spectrum in B —>■ X s , -f for which the 
leading distribution function drops out [§j||. In this case, even a model of the subleading 
distribution functions will provide a useful estimate of the theoretical uncertainty in these 
relations, and the resulting extraction of \V u b\. This investigation is in progress 0]. 
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